Breakdown of the cross-Kerr scheme for Photon Counting 
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We show, in the context of single photon detection, that an atomic three- level model for a transmon 
in a transmission line does not support the predictions of the nonlinear polarisability model known as 
the cross-Kerr effect. We show that the induced displacement of a probe in the presence or absence 
of a single photon in the signal field, cannot be resolved above the quantum noise in the probe. 
This strongly suggests that cross-Kerr media are not suitable for photon counting or related single 
photon applications. Our results are presented in the context of a transmon in a one dimensional 
microwave waveguide, but the conclusions also apply to optical systems. 
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The cross-Kerr effect, whereby the phase of one field is 
changed proportional to the intensity of another, arises 
from the nonlinear response of an atomic medium to ap- 
plied fields. It is usually described phenomenologically 
in terms of a third order term in the nonlinear polaris- 
ability, a description that is valid when the applied fields 
are strong and absorption is weak A derivation of 
the nonlinear polarisability description of the Kerr effect 
based on an underlying microscopic atomic model has 
been given by many authors including 0411] ■ 

Many proposed applications of the cross-Kerr effect 
however suppose that at least one of the fields is very 
weak — perhaps only a singlephoton — including non- 
demolition measurements quantum state prepa- 
ration [1-Qj]], quantum teleportation [l2| and quantum 
logic gates build-up (l3l - fl7l |. All these schemes require 
strong Kerr nonlinearities at the level of a single pho- 
ton. It is not clear that the standard model of a cross- 
Kerr effect, based on a third-order nonlinear polarisabil- 
ity, should be valid for fields with only a few photons. 

Doubts regarding the utility of the Kerr effect in single 
photon applications have been raised before. Shapiro and 
Razevi [H, took the multimode nature of the single 
photon pulse into consideration and found that there is 
extra phase noise compared to simple single mode cal- 
culations, leading to constraints on the achievable phase 
shifts. Gea-Banacloche [2fJ pointed out that it is impos- 
sible to obtain large phase shifts via the Kerr effect with 
single photon wave-packets. None of this prior work has 
addressed in detail the question of the cross-Kerr phase 
shift on a coherent probe field in the presence or absence 
of a single photon in the control field. 

Recently, superconducting circuits have become im- 
portant test-beds for microwave quantum optics, demon- 
strating quantised fields, artificial "atoms" (i.e. with well- 
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resolved energy levels), and strong "atom" -field interac- 
tions. The transmon [2l| is one of most promising super- 
conducting artificial atoms due to its insensitivity to 1// 
noise, strong anharmonicity, and large dipole moment. 
Indeed, the typical size of the transmon is comparable 
to the dielectric gap in an on-chip microwave waveguide, 
and so the dipole moment is within an order of magni- 
tude of the maximum that it can possibly be, given the 
geometrical constraints of the dielectric gap 22]. This 
fact leads to the possibility of very large cross-Kerr non- 
linearities, where the transmon provides the non-linear 
polarisability. Recent experiments using a superconduct- 
ing transmon in a ID microwave transmission line have 
demonstrated gigantic cross-Kerr nonlinearities: a con- 
trol field with on average 1 photon induces a phase shift 
in the probe field of 11 degrees [23| . Importantly, in this 
experiment, the microwave fields were freely propagating; 
no cavity was involved. 

This large cross-Kerr phase shift immediately suggests 
the possibility of constructing a broadband, number- 
resolving, microwave-photon counter, as long as the 
cross-Kerr induced displacement of the probe exceeds the 
intrinsic quantum noise in the probe. Indeed, broadband 
microwave photon counting is a crucial missing piece of 
the experimental quantum microwave toolbox, although 
there are several proposals for detecting microwave pho- 
tons dHH. 

In fact, the cross-Kerr interaction is strictly an effec- 
tive interaction based on weak field-dipole coupling ap- 
proximations. Ultimately it is mediated by the strong 
nonlinearities inherent in an anharmonic oscillator (e.g. 
an atom), so it must eventually break down as a useful 
description of the physics. It is therefore important to 
understand the contribution of the transmon (or atomic) 
dynamics to the effective nonlinearity in the limit of very 
strong coupling, which was achieved in (23j . In this work 
we investigate the coupled field-transmon dynamics in 
this limit, using proposals for microwave-photon count- 
ing as a technical objective to evaluate the validity of the 
cross-Kerr approximation. 
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FIG. 1: (Color online) (a) Illustrative experimental arrangement. 
A photon source emits a Fock state microwave photon into a ID 
planar transmission line with a S-type three-level transmon em- 
bedded in it. (b) Transmon level structure. The upper transition 
of the transmon is coupled with a coherent microwave field (probe) 
and its lower transition is coupled to the Fock state. The interac- 
tion induced phase shift in the probe field is detected by homodyne 
detection, (c) Cartoon of the Kerr-induced probe displacement. 



We consider two fields, a probe and a control, incident 
on a superconducting transmon qubit, which is treated as 
a three-level, S-type system in a one-dimensional trans- 
mission line. Such three-level systems are prototypes 
for analysing cross-Kerr nonlinearities We do not 
eliminate the transmon, but instead treat its dynamics 
exactly, including quantum noise in the incident fields. 
The probe is assumed to be a coherent field (or possi- 
bly squeezed), while the control field is in a Fock state, 
whose photon number, n, we are trying to measure. For 
our purposes, we restrict to n = or 1. 

We show that in spite of the very large cross-Kerr non- 
linearity, the induced probe displacement (i.e. the signal) 
in the presence of a single control photon is limited by 
saturation effects in the transmon, and is always less than 
the probe's own quantum noise. That is, the signal-to- 
noise ratio (SNR) is always below unity. Moreover, our 
conclusion also extends to the N-type four-level atomic 
level configuration, with which cross-Kerr media are of- 
ten modelled [29l - i32T | . These conclusions have rather pro- 
found implications for the exploitation of cross-Kerr phe- 
nomena in quantum technologies. 

The transmon levels are {\a), \ b), |c)}, with correspond- 
ing energy levels, taj, and decay rates, ji, as shown in 
Fig. [T] Relaxation between transmon energy levels is rel- 
atively fast compared to dephasing rates, which we ne- 
glect. The probe field, b, is in a coherent state |/3), and is 
nearly resonant with the \b) ■<-> |c) transition, whilst the 
control field is in a Fock state of n = or 1 photons, at 
a frequency u) C on close to the \a) <-> \b) transition. Qual- 
itatively, the control field induces a transient population 
transfer into the state \b), and the probe field induces a 
coherence, at, c , between states \b) and |c). This polari- 
sation couples back to the probe field, so that the probe 
field is modified from its input state according to the 



standard input-output relation 

bout = b ln + ^/Yc^bc- (1) 

The homodyne detector monitoring the output probe 
field yields a photocurrent given by 



(2) 



where y — —i- s / ? f^(a'bc — &cb) is the transmon polarisation, 
&ij = I*) 01 an d £dt = dW(t) is a Weiner process satisfy- 
ing E[dW] = 0, E[d 2 W] = dt. Finally, the useful signal 
is the integral of the homodyne current over the lifetime, 
T, of the photon wave packet 



fin 



dtJ%° m (t) 



(3) 



If n = the transmon dynamics are trivial, and £7[£o] = 
0. For n = 1, E[Si] ^ 0, and so Si represents the useful 
signal associated with a single photon in the control field. 
However, in any given measurement, the homodyne cur- 
rent includes quantum noise, characterised by the vari- 
ance (o\s,J 2 = Els 2 ] — E[S n ] 2 . To a good approximation, 
as n is independent of the photon number, n, and so we 
define the signal-to-noise ratio, SNR = E[Si]/ (y/2crs). 
Note that we assume that the homodyne current will 
also include technical noise sources. We ignore these, so 
that SNR represents the quantum limit for the proposed 
scheme. 

To study quantitatively the system consisting of a 
transmon interacting with propagating microwave fields, 
we adopt two different (but consistent) formulations, 
yielding both numerical and analytic results. 

In the first formulation we suppose the control pho- 
ton is generated by a fictitious cavity which is initially 
in a Fock state. The field in the cavity decays into the 
ID waveguide, and propagates to the transmon, which 
mediates the interaction between the control and the 
probe. We emphasize that the cavity is included sim- 
ply as a model photon source; the transmon is not con- 
tained within the cavity. To analyse this system, we use 
the stochastic cascaded master equation method [H, [34| . 
The stochastic master equation describing the condi- 
tional dynamics of the cascaded cavity field-transmon 
density matrix, p, is given by 



dp = (-i[H s ,p] + 7 con £> [a C on]p + D[Lb]p + T>[L c ]p)dt 
+ s/jcon([Lb, paly,,] + [a con p, L> b ])dt 

-™ /2 }pdW (4) 



H[L C 



where Lb = y/^bVab, L c = ^f c Obc and 

H s = A C (T CC + Ab&bb + Qp{<Jbc + <Jcb), (5) 

T>[f]p = — (2f pr* — pr* f — f p) , (6) 



H[r]p = rp + pf [ — Ti[ip + pi }p, 



(7) 



A h = uj ba - w con , A c = A p + A b (A p = uJbc - Up) and 
f2p = A /7 con /3. We solve Eq. @ for the conditional state 
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of the field-transmon system, from which we compute the 
conditional homodyne photocurrent, using Eq. @. This 
approach allows us to generate a simulated measurement 
record for ensembles of events in which n = or 1. From 
these simulated measurement records, we obtain a his- 
togram of homodyne currents, from which we estimate 
the SNR. 

The second formulation uses the Fock state master 
equation [H, Hfjj]. Instead of simulating the free space 
photon as the output of a fictitious cavity, the propagat- 
ing photon wave packet drives the transmon directly. The 
transmon density matrix acquires indices m, n represent- 
ing coherences between the transmon and photon Fock 
subspaces m and n. The corresponding master equation 
for the hierarchy of transmon density matrices p m ,n is 

Pm.n 

(t) = -i[H s ] + T>[L b ]p mtn + V[L C ] 

+Vnf*(t)[Lb,Pm,n-l] + Vmf(t)[pm-l,n, Ll\ 

where f(t) is a complex valued probability amplitude 
that determines the photon counting rate as proportional 
to |/(i)| 2 - We first solve the dynamics for po,o(i)j which 
drives Po } i{t) and pi${t), which in turn drives Pi t i(t). 
We solve these analytically and numerically, and use the 
quantum regression theorem [37[ to calculate the SNR 
(see part A of the supplementary information). The great 
advantage of this approach is that we obtain the SNR 
without resorting to stochastic simulations. 

If the photon is derived from the exponential (E) de- 
cay of a cavity mode, then f(t) = y/j con exp(-7 con i/2). 
Further, this method can handle arbitrary photon wave 
packets, and we include Gaussian (G) and rectangular 
(R) , as shown in Fig. [2jtop) . All pulses contains exactly 
one photon, that is J Q \f(t)\ 2 dt — 1 and their common 
width is l/7c „. The photon induces a polarisation in 
the transmon, shown in Fig. [2jbottom). We see differ- 
ent pulse shapes yield modest differences in the transmon 
polarisation (y{t)). 

In Fig. [31 the SNR is shown as a function of the probe 
amplitude with other parameters (the detunings and the 
single photon pulse width) optimised. The points are cal- 
culated by averaging 5000 trajectories of the stochastic 
master equation, whilst the solid line is computed from 
the Fock state master equation. There is good agree- 
ment between the two approaches. The inset shows a 
histogram of the results of the stochastic simulation at 
the value of /3 that optimises the SNR. Clearly the SNR 
is everywhere less than unity, and so it is impossible to 
reliably distinguish between zero and one photon in a sin- 
gle shot. The histogram confirms that the distribution of 
integrated homodyne current is much broader than the 
separation of the means. 

As an example of the influence of the various param- 
eters on the SNR, we plot the SNR as a function of the 
detunings At, and A c , as illustrated in Fig. [4] Clearly, 
the optimal SNR is located at Ab = A c = 0. 

The reason that SNR< 1 can be understood in the fol- 
lowing way: a single control photon induces a variation in 




FIG. 2: (Color online) The transmon responses for different 
control field wave packets, (top) blue dot-dashed curve: Gaus- 
sian pulse (G); Green solid curve: rectangular pulse (R); Orange 
dashed curve: Exponentially-decayed pulse (E); and (bottom) the 
corresponding polarisation response of the transmon. The pa- 
rameters are: A c = A;, = 0,-y CO n = 0.66727t,,7 c = 2^ b ,/3 = 
0.4 7 6(B)i 0.47 7 i,(i?); 0.59 7i) (G). 



the transmon polarisation b\, c , which manifests as a fluc- 
tuation in the homodyne current according to Eq. (|2|). 
However the polarisation of the transmon is a bounded 
operator: ||y|| < y/j b - The optimal photon wave packet 
width is T ~ 7^ (any shorter and the transmon can- 
not respond to the field; any longer and vacuum noise in 
the homodyne signal grows), so the integrated polarisa- 



tion can be no larger than < / Q T dt \\y\\ < j b ^ 2 . 

Quantum noise in Eq. ([2]) gives tr| > v&t[Jq dt£] = j^ 1 
we see that the signal to noise ratio is necessarily less 
than unity. Fig. [3] bears out this analysis: for small 
probe field amplitudes, the SNR increases quickly, how- 
ever the transmon dynamics quickly saturates at higher 
probe amplitudes. 

This argument suggests that the fundamental problem 
is the saturation of the transmon transition. It may be 
thought that this can addressed by increasing the number 
of transmons. We therefore briefly consider a system of 
N transmons, arranged such that the spacing between 
adjacent transmons is much smaller than the wavelength, 
the transmons are described by the collective atomic spin 
operators 



(9) 



The stochastic master equation describing the n- 
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H s = N{A C S CC + A b S b ) + ^/WcKSbc + S cb ) (11) 
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FIG. 3: (Color online) The SNR as a function of the probe ampli- 
tude j3 at optimal parameter setting. The orange square represents 
the numerical SNR from the stochastic, cascaded master equation 
method and the green curve represents the analytical SNR from 
the Fock state master equation method and quantum regression 
theorem. The inset provides the histogram of the highest SNR for 
zero and one signal photon. 



SNR in the A -A c plane 




FIG. 4: (Color online) The SNR of the probe field as functions 
of detunings A;, and A c . The other parameters are: 7 C = 2 7 j,, 
7co „ = 0.6772 7i) , = 0.47;, . 



transition system is given by 

dp = -i[H s ,p]dt + j con V[a con ]pdt + N%T>[S ab ]pdt 

+ Nj c T>[S bc ]pdt - V Nj con j b ([S ba , a c p] + \pa\, S ab ])dt 
+ ^fW c H[S bc e-™' 2 ]pdW (10) 



We can see that the ensemble master equation (|T0|) is 
the same form as the single-transmon master equation, 
albeit with decay rates and energies scaled by a factor of 
N , and its dynamics must therefore be correspondingly 
faster. However, this cannot change the optimised SNR, 
so the SNR for the N transmon case will be the same as 
for the single transmon case. 

It is worth commenting on a number of other avenues 
that we have explored, but which yield similar negative 
results. Firstly, it may be thought that squeezing the 
probe field in an appropriate quadrature would reduce 
the noise, and therefore improve the SNR. Since we are 
monitoring the displacement of the probe field, we should 
squeeze the phase quadrature. However this will enhance 
noise in the conjugate, amplitude quadrature. The ad- 
ditional noise in the probe amplitude adds noise to the 
transmon dynamics arising from fluctuations in Q p . We 
find numerically that these tradeoffs yield no net im- 
provement in the SNR. (see part B of the supplementary 
information) 

Secondly, there is another multi-transmon limit, in 
which the transmons are sufficiently separated that they 
may be considered as a series of cascaded systems. It 
might be hoped that the control field will interact sequen- 
tially with each transmon. A local probe at each trans- 
mon would then yield an independent estimate of the 
probe displacement. After M such independent probes, 
the SNR would be improved by a factor of M 1 ' 2 , and 
the SNR could be made arbitrarily large by increasing 
M. However, the Kramers-Kronig relations impose a 
tradeoff between the phase shift at each transmon and 
the probability that the control field is reflected: a large 
phase shift necessarily implies a large reflection probabil- 
ity. Again, we find numerically that the tradeoff yields 
no net improvement (see part C of the supplementary 
information) . 

Thirdly, a number of proposals for inducing cross-Kerr 
nonlinearities in optical systems use an N-type four-level 
system [29|, l30J , with a strong classical field addressing 
the intermediate transition. As we show in part D of the 
supplementary information, in the limit of strong driving, 
this maps onto the same three-level structure we consider 
in this work, so the conclusions we have reached here also 
apply to such N-type systems. 

Fourthly, We also numerically investigated the effect 
of varying the ratio ~f c /lb arid for 1 < "f c /jb < 100 find 
that the SNR is still much less than unity (see part E of 
the supplementary information). 

A number of proposals suggest using weak Kerr media 
to build controlled phase and C-NOT gates with fewer 
resources than linear optical schemes [3, [ljj ■ In these 
schemes the cross-Kerr phase shift per photon is much 
less than n, so a strong coherent bus compensates for the 
weak nonlinearity, such that the small cross-Kerr phase 
shift manifests as a large displacement of the strong co- 
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herent field. However the saturation of the cross-Kerr ef- 
fect described above indicates that once the displacement 
of the strong coherent field approaches its own quantum 
noise, saturation effects lead to the breakdown of the ef- 
fective cross-Kerr description, rendering such protocols 
ineffective. 

In summary, we have investigated the feasibility of mi- 
crowave photon-counting based on an induced cross-Kerr 
nonlinearity arising from coupling to a large anharmonic 
dipole. We find that saturation of the transmon transi- 
tion limits the SNR to less than unity. As such, it is not 
possible to use strong, atom-induced cross-Kerr nonlin- 
earitics to perform single photon detection. This conclu- 
sion applies to a number of extensions of the basic model, 
including multiple transmons, cascaded transmons and 
N-type, four-level system. Further, it limits the applica- 
bility of any proposal that requires a cross-Kerr nonlin- 
earity to produce a displacement of a coherent field by 
an amount greater than the intrinsic quantum noise in 
the coherent field: it is precisely this condition where the 
effective cross-Kerr description breaks down, and satura- 
tion effects become dominant. 

This work is sponsored by ARC Center of Excellence 
for Engineered Quantum Systems and the China Scholar- 
ship council. Bixuan Fan would like to thank Dr.Zhenglu 
Duan and Dr. Matt Woolley for helpful discussions. 



SUPPLEMENTARY INFORMATION 

A. The analytical solution for a three-level system 
on resonance 

Here we provide the details of the analytical solution 
for a three level system. As presented in the main text, 
we consider a three-level transmon coupling with a coher- 
ent field at the \b) •<->• |c) transition and a single photon 
at the | a) o |6) transition. The one-photon Fock state 
master equations (35j are given by 

A>,o(*) = -i[H s , POt o]+V[L b ]p 0t0 + T>[L c ]pofi (12) 
A),i(t) = -i[H s ,po,i}+V[L b }po,i+V[L c }p s (13) 

+ f*(t)[L b ,Po,o) 
pi,i(f) = -i[H s ,p ltl }+V[L b ]pi tl + V[L c ]p lt i (14) 
+ r(t)[L blPl , Q }+f(t)[p 0A ,Ll} 

where the temporal profile function f(t), system opera- 
tors L b , L c and the system Hamiltonian H s have been 
defined in the main text. Here p m ,n is a matrix and the 
subscripts m and n denote the photon number basis. 

Initially the transmon is prepared at the ground state. 
The lowest equation for po,o(f) can be easily solved as 
Po,o(t) — Po,o(0). Then it is substituted to the next equa- 
tion for po,i(i)i which is traceless. For our system and 
an arbitrary input Fock state, the generalized density 
matrices p m ,n in the Bloch-like representation can be pa- 



rameterized as 



-16(m, n) + -a„ 



(15) 



where I is a 3 by 3 identity matrix and vectors a mj „ = 

(^lmni ^2mni ^3mni ^4mnj ^5mni ^6mnj &7mri) Og mn ) and 

A = (Ai, A 2 , A3 , A4 , A5 , Ag , A7 , A 8 ) . The Gell-Mann 
matrices [381 ] for qutrit are 



Ai 



A 4 = 




A 2 = 



'0 1 

I A 5 = 

1 





A 3 = 



'0 CP 
A R = I 1 

,0 1 0; 



'0 
A 7 = ( -i 
,0 i 




(16) 



These matrices are traceless, Hermitian, and obey the 
relation tr(AiAj) = 2§ij. 

Substituting po,i m terms of ao,i into the Fock mas- 
ter equation, we can have the coefficient equations for 
a>iOi(t)(i = 1,2,. ..8). With initial conditions Oioi(0)(i = 
1, 2. ..8) = 0, we have the solutions: 

0101(f) = a 20 i(t) = 0301(f) = o 8 oi(f) = (17) 

0501(f) = CiC 5 exp(-0 1 f)+CiC 6 exp(-0 1 f) (18) 

+ C*iC 7 exp(-7 co „f/2) 

o 60 i(f) = CiC 2 exp(-0 1 f)+CiC 3 exp(-0 2 f) (19) 

+ CiC 4 exp(- 7co „f/2)] 

a 4 oi(f) = «O50i(f) (20) 

0701(f) = -io 6 oi(f) (21) 



3 7b /4 + 2 = 3 7b /4 - ^0/4 

V-32/3 2 + 7b 



C 2 

c 3 

C 5 
C 6 

cv 



\fib~Tc 



where 

(22) 
(23) 
(24) 

0[2 7b (4/32+ 7b )-3 7b7co „+ 7c 2 o J V ' 

4V7b"(8/? 2 -7b + VTbO) + lcon {^T b - 9) (25) 
4 % / 7b ~(-8/3 2 + 7b + ^T b 6 - 7co „(V76 + 9) (26) 
-4 lb 9 + 2 lcon 9 (27) 
2^2/3(-3 7b + 2 7co „ + ^f b 9) (28) 
272/3(3 7b - 2 7con + ^f b 9) (29) 
4 v /2 7 7/30 (30) 

and the density matrix po,i at time f can be presented as 



Po,i(f) = 77 a O,l A 







a 40 i(f) a 60 i(f) 



(31) 
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where we used 7 C = 27^ (for a transmon). 

The top level equation pi t i(t) which represents the ac- 
tual system evolution is not traceless. In the Bloch rep- 
resentation the equation is pn = |l+ \a\ ,iA with initial 

conditions: a%n(t = 0) = — 2/\/3 and am(t = 0) = for 
i ^ 8. Substituting this expression and the solution for 
Po,i hito the master equation we have the motion equa- 
tion for ai i. 

Only motion equations for 0211, 0311 and asu are cou- 
pled and non-zero. We define a vector: 



x = (a2u,a3ii,a8ii)' 3 



with initial condition 



z(0) = (0,0,--=) T 



Tt =A X + B{t) 



with coefficient matrices: 

-3 7 fc/2 -2^/3 
.1=12^277/3 -5 76 /2 -^76/2 
y3 76 /2 - 7b /2 



(32) 



(33) 



(34) 



(35) 



~2\/7ft£a50i 

S(t) =| -76 + 2^/6^601 
-^-2^376^601 



(36) 



By diagonalization and integration, we obtain the solu- 
tions as: 

x[i]{t) (37) 
= VnC xl {t)[ J [Cni-t'WnB^f) + C n (-t')Q 12 B 2 (t') 

+ Cu(-t')Qi 3 B 3 (t')]dt' - ^=Q 13 ] 

+ V 12 C 22 (t)[ f [ai2(-i')Q 2 iBi(0 + C 22 (-t')Q 22 B 2 (t') 



where 

Si = Cii(i) JiCui-t^QMt^ + Cui-t^QMt') (38) 
+ C u {-t')Q i3 B 3 {t')]dt' 
= ^ ^Q i2 + -^=Qi 3 ^ (1 - cxp (AJ)) 

^JlblconC\Ci r- 

+ —t (Qa - v3(y i3 )(exp(-(7b + 7 co „)i) 

7con T" 

- exp(Ait)) 
2^767 

C1C2 



J-J (0 i2 - \/3Q i3 ) (exp(-(A! 



Ai + Ai + — 

- exp(Aji)) 

- (g a - V3Qis) (exp(-(A 2 

^2 + \ + 2 ~ K ' 

- exp(Ajt)) 



0<) 



with matrices V , Q and C being the eigen-vector matrix, 
inverse eigen-vector matrix and the diagonalized matrix 
of the coefficient matrix A. 

Then we have the system density matrix at time t: 

1, 



(39) 

\ 


qgii(*) 

/ 

The transmon polarisation for the homodyne detection 



(a 311 (t) + ^f -ia 211 (t) 

ia2ii(t) " ss ^-a sll {t) 
y 



+ c 22 (-t')Q 23 B 3 (t')]dt' - - 7 =g 23 ] 

+ V 13 C 33 (t)[ /"[Ci3(-t')Q 3 iBi(0 + C 33 (-t')Q 32 B 2 (t') 



C 33 (-t')Q 33 B 3 (t')]dt' - -j=Q 33 ] 

V a Sx + V l2 S 2 + V i3 S 3 - -^VaCt^Qu 

-^V l2 C 22 {t)Q 23 - l=V i3 C 33 (t)Q 33 



(y(t)) = Tr[-i(L c - L,t)p lfl ](t) = Tr [A 2 pi,i] (t) = a 211 (t) 

(40) 

In Fig. [5l we compare the results for (y(t)) by analyt- 
ical and numerical methods. It is a perfect agreement 
between the results by the two methods. 

The noise or the variance of the detected signal, can 
be calculated by the quantum regression theorem [37[ : 

(AS) 2 = E [S 2 ~ S 2 ] (41) 

dt [ di'u(i'-i)Tr[y(t)e £(t '- t) (-iL c p(t) 
Jo 

+ip(t)Lt)] + u(t - OTrtye^^-iLcKt') 
+ip{t')L\)]+T-S 2 

where 5 = fT dtJ% om (£), S — J? dty uc (t) and the func- 
tion u(t) = l(t > 0);u(t) = 0(i < 0). The subscripts uc 
means unconditional results. 

B. Squeezed probe 

In this subsection we replace the coherent probe field of 
a phase-squeezed state and the corresponding stochastic 
master equation is 
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FIG. 5: (Color online) The time evolution of the transmon 
polarisation: comparison between analytical and numerical 
results. The parameters are: ~/ c = 2^ b ,^ con — j b , A c — A b — 
0,P = 7b. 



FIG. 6: The SNR as a function of the squeezing degree (db). The 
parameters are: 7 C = 2-y b ,^ con = 0.6672~f b ,/3 = 0.47 6 , A c = A 6 = 
0. 



dp = (-i[H s ,p] + ^ con V[d con ]p + j b T>[a ab ]p 

\Jlconlb{[&ba, 0,conP\ + [pO-con' &ab 

"f c NV[a cb ] + "f c Ma bc pa bc - -f c M*a cb pfr cb )dt 

-m/2 _ ^ N 



l)V[a bc ] 



low dimensional systems photons will be reflected com- 
pletely by qubits on resonance [IlJjIciJ. On the other 
/.rjxhand, EIT-like effects appear when there are more energy 
levels [ijj ■ In this three-level structure, two-pathway in- 



^U[(N + 1 + M)a bc e 



M*)a cb e 



where 



H s = A c a cc + A b a bb + ^pf c j dvf3 v (& bc + & cb ) (43) 

with M = sinh(r) cosh(r)e* 6 ', N = sinh 2 (r) and L = 
1 + 2 N + M + M* . The instantaneous photocurrent is 

I h c ° m it) = (y) c (t) + VL£(t) (44) 



terference forms and the familiar EIT-like transparency 
window appears, shown in Fig. [7J However, how wide is 
the window and whether the parameters for this trans- 
^JpdWparency window is consistent with those at best SNR are 
questions. In the following we calculate the transmission 
rate by the similar procedure in [39j . 

The cascaded system is equivalent to the model of 
transmons interacting with a single photon pulse cen- 
tered at uj c and a coherent field in the transmission line 
directly. After linearizing the dispersion in the vicinity 
of uj c and formally adding the non-Hermitian damping 
terms (which come from the Markovian approximation 
by tracing out the bath operators), the Hamiltonian in 
the real space can be written as 



This equation indicates that the noise term £ is multiplied 
by a factor of \f~L. Notice that when 9 is chosen as zero, 
L = e~ 2r and the noise is reduced while when 6 = ir, the 
noise is amplified. 

In Fig. [S] we show that the squeezing in the phase 
quadrature can only help to improve the SNR slightly. 
According to uncertainty relation the phase squeezing 
indicates an amplification of the amplitude noise, which 
results in larger dynamical noise in the atomic response. 



H 



dx[d} R {x)(uj 

(uj b - ij b /2)a bb 
a' L a ab + o- ba a L ] 



d 



)dn(x) 



d_ 

9 di 



ug dx' 
)a L {x)} + {u> c 



(45) 

,, - ij c /2)a cc 



dx^f%5{x)[a , R (T ab + a ba a R 
~ c a{<7 bc 



CTcb 



C. Cascaded n transmons 

In the main text we have discussed the single transmon 
and a transmon ensemble for single microwave detection. 
In this subsection, we investigate the feasibility of cas- 
cading multiple transmons, each with a probe field. 

First we evaluate the transmission rate of the signal 
field after passing one transmon. As is known that in 



where v g is the group velocity of the signal photon, which 
depends on the geometry and material of the waveguide. 
For a typical coplanar waveguide of quantum circuit sys- 
tem, v g = ^j^, L , — cF% = c l \ /6e ff w ith the effec- 
tive permittivity e e ff around 5.9. The time-independent 
eigen-equation is 

H\E k )=E k \E k ) (46) 
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with Ek = u) = lj c + v g k R and 

\E k ) = Jdx(f> R a R [x) |0fl,0i,o) (47) 

+ J dxc/) L d} L \0 R ,Q L ,a) 
+ a \0 R ,0 L ,b) +c 2 \Qr,O l ,c) 
Then we have the equations for the coefficients: 
d 

(ujcon ~ WbI-q^&R + 9\5{x)ci = lo4> r (48) 



(uj con + iv g)-Q^ ( t ) L + ad(x)ci = u<j> L (49) 




-5 -4 -3 -2 -1 1 2 3 4 5 



(u) b - *y)ci + y/jb(^ R (0) + 0l(O)) + V> c 2" = wc{50) 



where 



(uj c - - w P )c 2 + y/^cia — loc 2 (51) 

<j) R = exp(ikx)8{—x) + texp(ikx)6(x) (52) 

</>l = rexp{— ikx)8{— x) (53) 

with 9(x) being the Heaviside step function. 

By solving the equations above, the transmission am- 
plitude can be obtained; 



FIG. 7: (Color online) Transparency spectrum. 

Probe 1 Probe 2 Probe n 



Signal 



4 



FIG. 8: (Color online) The schematic for microwave photon 
counting using cascaded n tranmons. 



t 



(A' h +i^r c /2){A' c + i lc /2)- lc o? 
(A' b + i lc /2 + ij b /vg)(A' c + i lc /2) - 7c a 2 



(54) 



where A' b = uj b — w and A' c = uj c — lj p — w. Here the fre- 
quency width of the signal photon is much smaller than 
its central frequency, that is to say, it is a relatively nar- 
row pulse, therefore in the the following we use A& and 
A c to replace A' b and A' c . 

From Fig. an induced transparency window appears 
in our system as expected, which is caused by the two- 
channel interference between transmon transitions. The 
width of window is twice of the coupling ,J^%a. One can 
find that to achieve a high transmission rate, either large 
a or large Ab is required. This indicates that a large 
transmission rate corresponds a low SNR. 

Assuming we cascade n transmons with separate 
probes and detectors, seen in Fig. [8] ideally it is equiv- 
alent to average over n trajectories. However, after in- 
cluding the reflection, the effective SNR n becomes (high 
order terms o(R) are omitted here): 

n-1 

SNRn = SNR 1 (y/nT n ~ 1 + ^2j/V^T j ~ 1 R) (55) 

3=1 

where SNRi is the SNR in one transmon and one probe 
case. After optimizing parameters, we obtain the SNR 
as a function of probe number as shown in Fig. [SJ Obvi- 
ously, the signal can not win the noise in this method. 



D. Conversion from N-type Four- level structure to 
ladder-type three-level structure 

The N-type four-level structure has been suggested to 
be a promising candidate for implementing cross-Kerr 
nonlincarity. Here in this subsection, we show that it can 
be approximately mapped to a three-level ladder system. 




5 10 15 20 25 

Number of cascaded tranmons n 



FIG. 9: The plot of SNR as a function of transmon number 
u. 
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I 3 ) 



P 



FIG. 10: Illustration of the conversion from an N-type four- 
level system to a ladder-type three-level system. 



The Hamiltonian for a four-level system coupling with 
a signal field ft at transition 0-1, a control field ft at 
transition 1-2 and a probe a at transition 2-3 in a rotating 
frame is given by 

H = (Aio - Ai 2 )a- o + A120-11 + A32633 (56) 

CI, 



+ \/7oi/9(a'io + <7oi) + -^-(012 + 021) 

+ V7 32 Q; ( (3 '23 + CT32) 

When the transition 1-2 is strongly driven by f2, the sub- 
system of the transition 1-2 and the control field can be 
diagonalized independently. We define a unitary trans- 
formation U as 



cost 
sin( 



— suit 
cos (9 



By diagonalizing U^Hi 2 ,nU, we have: 
9 = — arctan(51/Ai2) 



(57) 



(58) 



and the states |1) ,12) can be represented in the dressed 
states I—) and |+)|42|: 



|1) = cos0|-) -sin0|+) 
|2> = cos0|+)+sin0|-) 

Then the Hamiltonian can be rewritten as 



(59) 



As shown in Fig. [IT] it is a good approximation when 
there is a strong driving (57 3> 701, a, ft) at the dressed 
transition. 



— N-type four level 

- Effective ladder three level 




FIG. 11: (Color online) The comparison of the population of 
the excited state in the two configurations. The parameters 
are A 12 = 0, A 10 = A 32 = A_, fi = IO710, and a = ft = 
712 = 732 = 710 • 




FIG. 12: (Color online) The SNR as a function of the ratio 7c/7b- 



H = — Aioo-qo + A32033 + A+(7 ++ + A_<7 — (60) 
+ V7oi/3(cos(9((T_o + 5-q-) - sin0(<7 +o + a +)) 
+ % /7 32 a(cos(5( l T + 3 + 0-3+) + sinf5(,T_3 + <t 3 _)) 



where A± = |A 12 ± \y/ + 

We assume the signal field and probe field are tuned 
to be resonant with the state |— ), that is, Aio = A_ 
and A 32 = A_. We perform a rotating and ignore the 

fast- varying terms like eW A ~ 12+n ' t , then we have an ap- 
proximate three level system: 



H 



/7oT/3((T_o + o-q-) + sin(9 v ^yi2"a((T_3 + 0-3 -061) 



E. Effect of different relaxation rate ratio "fc/'ft 



For a transmon, the relaxation rate of the upper tran- 
sition is twice that of the lower transition but in other 
three-level system we may have different ratio of the re- 
laxation rates. Here we investigate the effect of different 
relative relaxation rates on the SNR numerically. Fig. 
[T2l shows the SNR as a function of the ratio jc/jb- It is 
clear that no matter lower ratio or higher ratio the best 
SNR can not be improved much and still far below unity. 
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